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ABSTRACT: It is expected that the only effect of heavy dynamical fermions 
in QCD is to renormalize the gauge coupling. We derive a simple expression for 
the shift in the gauge coupling induced by N j flavors of heavy fermions. We 
compare this formula to the shift in the gauge coupling at which the confinement- 
deconfinement phase transition occurs (at fixed lattice size) from numerical simu- 
lations as a function of quark mass and Nt. We find remarkable agreement with 
our expression down to a fairly light quark mass. However, simulations with eight 
heavy flavors and two light flavors show that the eight flavors do more than just 
shift the gauge coupling. We observe confinement-deconfinement transitions at 
[3 = induced by a large number of heavy quarks. We comment on the rele- 
vance of our results to contemporary simulations of QCD which include dynamical 
fermions. 



1. Introduction 

QCD investigations frequently deal with the effect of heavy fermions either as real phys- 
ical effects (heavy quarks) or as the consequence of the regularization (Wilson fermions). 
In all cases the influence of heavy fermions at low energies was expected to be no more 
than some induced effective gauge coupling. 

Finite temperature simulations with two light and one heavier quarks do not show a 
significant difference from two light quark simulations. The effect of the heavy fermion 
can be described to a good approximation by a shift of the gauge coupling, A/3 ~ 0.08 for 
m q = 0.25 (Ref. 1) and A/3 « 0.13 for m q = 0.1 (Ref. 2). 

Wilson fermions have 15 heavy fermion doublers for each light fermion yet the spectrum 
hardly differs from the spectrum of staggered fermions if one takes into account the doublers 
by shifting the gauge coupling. For two flavors at a gauge coupling around (3 = 5.6 and 
hopping parameter value k « 0.16 this shift is about A/3 « 0.3 - the only apparent effect 
of the doublers. 

When can we expect that the fermions influence the physical spectrum in a non-trivial 
way and when can we just replace them with an effective local gauge action? The answer 
obviously depends on the physical processes we are investigating. Heavy fermions are 
always present in the spectrum, unless their mass is above the cut-off, but if the low lying 
gauge and light quark hadronic spectrum is much below the energy level of the heavy 
fermions they will not directly influence the low energy spectrum. 

The fermions' induced gauge coupling can be calculated by evaluating a 1-loop graph 
if the fermions are heavy. This analysis was presented in Ref. 3 using dimensional regular- 
ization, where the possibility of generating a continuum gauge theory with heavy fermions 
was investigated. The lattice regularized calculation is briefly mentioned in Ref. 4. In this 
paper we analyze further the analytically predicted induced gauge coupling and compare 
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it to existing and new numerical results. 

2. The induced gauge coupling on the lattice 

Consider the lattice regularized model of iVj fundamental (Wilson) fermions interacting 
with SU(3) gauge fields, whose action is 

S = (3 V Tr(U p ) + ±- V ^ n Knm[U]^ m , (1) 
n,/i n,m 

where 

^nm[^ = ^nm ~ K ^ ~ l/fi}Unfj,^n+iJ,,m ~H ( r ~t~ lf/j,)Unfi^n—fi,m)- (2) 
k is related to the inverse of the bare fermion mass 

« = " , (3) 

2ma + 8r v y 

where a is the dimensional lattice spacing, r = 1 corresponds to the usual Wilson fermion 
formulation while r = describes Nf = 16 x N f staggered fermions. Integrating out the 
fermions we obtain the effective gauge action 
S e// = S 5 -TrlnK[U] 

(4) 



= S a + 



where the sum is over all closed gauge loops V and l[T] is the length of the loop. Using 
the continuum representation of the gauge field U n ^ = e ia 9^A n ) one can express i n 
terms of the continuum fields A^(n) as the sum of one loop diagrams 



3 



The leading term of the effective action is the usual continuum gauge action -^F^yF^y 
where the coefficient 1/<7q can be calculated by evaluating the two 2-legged graphs in Eqn. 
5. The quantity 1/<7q can also be calculated starting with Eqn. 4 and using the method 
presented in Ref. 5 for adjoint scalars. It happens that this technique is actually incorrect 
for adjoint fields but correct for fundamental ones. 



The result is given by a four-dimensional lattice integral 
where S(p) is the lattice fermion propagator 



S = — -rJ^cos(p M ) -z^7 M sin(p M ) (7) 

and Q(p) is given by 



2k 



Q{Pn) = ir sin(p M ) + 7(U cos(p M ) . (8) 
The integral reduces to the hopping parameter expansion result in the k — > limit 

\ = AN f n A , r = 1, 

? (9) 

-s = 2N f K 4 , r = 0. 

#0 



For small ma (k — > 0.125) it has a logarithmic singularity 

1 Nf 7T 2 

The effective action has additional terms containing more derivatives and/or external gluon 
legs. These graphs are multiplied by negative powers of m and are suppressed for heavy 
fermions 3 . In the limit where the higher order terms can be neglected the effective action 
is indeed a pure gauge action with bare coupling constant given by Eqn. 6. In terms of the 
plaquette action lattice model it corresponds to an effective plaquette term with coefficient 
A(3 = 6/(7q. Table 1 shows A/3 for several m values for r = and Nf = 16iVy = 1 fermion 
flavor. 

3. Validity of the effective action 

In this section we investigate under what conditions a single plaquette effective action 
can describe the fermionic theory at low energies. 

The question is two-fold: 1) Can the non-local effective action Eqn. 4 indeed be 
replaced by a single plaquette term and 2) how well does Eqn. 6 predict the coefficient of 
this term? It is possible to have a pure gauge effective action in a region where Eqn. 6 is 
no longer valid. We will consider the second point in the next chapter and now investigate 
the first question. 

The low energy effective theory can be considered pure gauge if the gluonic spectrum 
characterized by the A parameter is much lower than the fermionic mass scale that is 
characterized by the fermion mass m. 

The one-loop definition of the lattice A parameter is 

A '-4 exp {-if}- < n > 

where [3q = HN c /48n 2 is the first (universal) coefficient of the (5- function and j3 e ff = 
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(3 + A/3. The condition A; att « m can be expressed as 

^y^ + ln(am)>>0. (12) 

As A/3 is proportional to Nj, this condition can be translated into a lower limit on the 
fermion flavors. For example, for (3 = 5.7, ma = 0.1, assuming the validity of Eqn. 6, if 
Nj » —28 (i.e. for any physical Nf) the heavy fermions and the gluonic sector decouple. 

It is interesting to consider the (3 = strong gauge coupling limit. For small m A(3 is 
logarithmically divergent leading to the condition for decoupling 

N f > | (13) 

in the a — > limit (Recall Nj = 16Nf). The minimum number of flavors coincides with 
the value where the (3 function of the gauge-fermion system changes sign and becomes non- 
asymptotically free. Since the presence of a gauge term {(3 ^ 0) relaxes the limit on Nf, if 
Eqn. 6 holds in an SU(3) gauge theory with 17 or more fermions, then the fermions always 
decouple from the low lying gluonic spectrum. An identical condition was found in Ref. 3 
using dimensional regularization. An interesting consequence is that the m = theory is 
always deconfined even in the strong gauge coupling limit for Nf > 17 flavors, assuming 
A/3 diverges as in Eqn. 10. One should keep in mind, however, that this derivation is valid 
only if the higher order terms in the effective action can be neglected. 

4. Examples 

Now we want to address the second question. When a gauge-fermion theory with 
massive fermions can be replaced by a gauge theory with a single plaquette action, is the 
shift in (3 induced by the fermions given by Eqn. 6? The easiest way to explore the shift in 
(3 due to fermions is by tracking the confinement- deconfinement transition as a function 
of quark mass and number of flavors. 

The quenched phase transition at Nj< = 4 is at (3® = 5.69(1) [6]. Introducing N j flavors 
of fermions with mass m will shift the transition to (3^ f = (3® — A(3. If m is such that 



the fermionic action can be considered pure gluonic at low energies then Ap(Nf,m) = 
NfA(3i(m). If, in addition, m is large enough that the perturbative formula is valid, 
Af3\{m) is given by Eqn. 6. Thus we expect the following behavior for the shift Af3(Nf, m): 
For m » A where the fermion and gluon mass scales are well separated we expect to 
see universal behavior A/3/Nj- = f(m) where f(m) is given by Eqn. 6. For smaller m 
we expect Eqn. 6 to fail quantitatively. However, it might happen that Af3/Nj is still 
some universal function of the quark mass. Finally, when the fermion scale is the same 
order as the gauge scale one can no longer replace the fermions by an effective gauge 
action. The shift Af3/Nj would then be different for different Nj, N T . Measuring the 
finite temperature transition for different Nj- and m values makes it possible to distinguish 
the different scenarios. 

The finite temperature transition is first order for the pure gauge theory, and is stable 
under the inclusion of heavy fermions. With decreasing quark mass the location of the 
transition shifts downward in (3. At some point the deconfinement transition line termi- 
nates (at sufficiently light quark mass) . We might still be able to track the crossover point 
as a function of Nj- and m. As long as the fermionic spectrum remains heavy compared to 
the low energy gluon spectrum, the system could still be described by an effective gauge 
action and Eqn. 6 could be valid. 

At very small or zero quark mass (depending on the number of light flavors) there is 
a second transition whose behavior is thought to be primarily chirally-restoring. At this 
transition the role of the fermions is fundamental and one would not expect the decoupling 
of the gluonic and fermionic spectrum. 

In a model with light flavors of mass mi and heavy flavors of mass m^, we also 
expect that the transition should be shifted by the heavy flavors: /? c (iVj,ra;, A^m/J = 
f3 c (Ni,mi) + AP(Nf l , mfj) where A/3 is given by Eqn. 6. It is a phenomenologically in- 
teresting question to ask, "How light is still heavy?" For example, several groups have 

7 



recently performed simulations with iVj = 2 and iV^ = 1 in an attempt to model the de- 
confinement transition in the real world of two light (u, d) quarks and one strange quark. 
To the extent that Eqn. 6 predicts the shift in lattice critical coupling, the heavy flavor is 
merely renormalizing the gauge coupling and contributing no new physics. 

In the above consideration we had to assume the relation f3 = 6/<7q - the induced gauge 
coupling is expressed through the bare continuum coupling g^ while in a lattice simulation 
one uses the coefficient of the plaquette term f3. (3 = 6/<7q should hold in the continuum, 
large (3 limit; one expects to encounter deviations when the finite temperature transition 
happens in the strong coupling (small (3) region. 

The A/3 values of Table 1 are calculated on an infinite lattice. On a finite lattice 
simulations there are corrections to it. For example, at iVy = 4 the fermions induce an 
explicit Polyakov loop term in the action which is the same order as the induced (3 of the 
gauge coupling (k 4 in hopping parameter expansion). In all simulations with dynamical 
fermions at iVy = 4, we have a nonzero expectation value for the Polyakov loop, even in 
the confined phase. That is an additional source of error in comparing Eqn. 6 to finite 
lattice numerical results. 

Now we consider a number of cases. We have chosen to focus on staggered fermions 
since it is easier to make a connection to Eqn. 6 with them than with Wilson fermions. In 
Sect. 5 we discuss a possible connection to the confinement-deconfinement transition for 
Wilson fermions. 

In the following we will translate numerical data to express the shift in the gauge 

coupling caused by one of the fermions only, Aft = (/$ - f3c f )/N f . Here $ is the Monte 

Nf 

Carlo quenched critical coupling and (3 C is the Monte Carlo Nj flavor critical coupling. 
This way we can compare simulations with different Nj and iVy values. 
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4.1 N f 



= 24 



Table II shows the result of a 24 flavor staggered fermion simulation on 6 d x 4 lattices 
for several mass values. These simulations were done by us and use a version of a code 
written by the MILC collaboration^ 7 . We employ the hybrid molecular dynamics algorithm 
described in Ref. 8. We have defined the dynamical fermion fields on all sites of the lattice, 
so that the "natural" number of flavors in the simulation is a multiple of 8. Simulations 
with large Nf require a very small timestep compared to ones with small Nf since the 
fermion force in the microcanonical evolution equations scales linearly in Nf. For example 
for m = 0.5 stepsize At = 0.020 is needed. The Nf = 24 transition is very sharp. Fig. 
1 shows the expectation value of the Polyakov loop at m = 0.5. At (3 = 4.62 the time 
evolution shows tunneling between two states (Fig. 2), the transition is probably first 
order. The data points agree with the analytical prediction for m > 0.25. At m = 0.25 
simulations with Nf < 8 agree with the analytic prediction. The deviation here should be 
attributed to the fact that j3 c = 3.90(5) is a very strong coupling where [5 = Q/g^ does not 
hold anymore. 

4.2 Nf = 17 

These data (shown in Table III) are from runs using the Langevin updating algorithm 
on Nt = 4 lattices? The analytic formula consistently overestimates the shift in (3 C . This 
is hard to understand given that the Nf = 24 and Nf = 8 simulations (see below) are 
well represented by the formula. However, the Langevin timestep Ati is related to the 
timestep of microcanonical simulations Atj^ by 10 At^ = (At M ) 2 /2. These simulations 
are performed at At^ = 0.01 corresponding to At^ = 0.14, which is known to be large 
enough to induce sizable integration timestep errors. 
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4.3 Nf 



= 8 



These data, shown in Table IV, are also from runs using the Langevin updating al- 
gorithm on iVy = 4 and 6 lattices. For smaller mass values the results are very sensitive 
to the step size used in the simulations. For too large At the transition is generally over- 
estimated, so the shift A/3 is underestimated. The analytic formula accurately predicts 
the location of the transition or crossover point for the larger values of the quark mass 
studied. For smaller quark mass values the agreement is still reasonable though Af3 M< ~ ! is 
consistently smaller than the analytic prediction. At very small m and iVy > 8 a number 
of authors 11 ' 12 ' 13 have seen a transition which may be a bulk transition. Our analytic 
formula does not predict this transition. On the other hand at such small mass values 
Eqn. 6 is not expected to be valid anymore. 

4.4 Nf = 4 

These simulations, shown in Table V, do not show a phase transition at moderate 
values of the quark mass. At small m they show a first order transition which is believed 
to be associated with chiral restoration. The location of the transition/crossover is well 
predicted by Eqn. 6 down to m = 0.05. For ATy = 4 at m = 0.073 the first order chiral 
transition switches on? It is surprising that Eqn. 6 is still valid. One can see deviation 
from the analytic formula for m < 0.025. 

4.5 Nf = 2 

Most of the Nf = 2 simulations were performed at very light values of the quark mass. 
They do not show a phase transition; instead, they show a smooth crossover from a chirally 
broken phase to a chirally restored one. Nevertheless, the location of the crossover point 
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is very well tracked by the analytic formula, even at very light values of the quark mass. 
The results are collected in Table VI. 

4.6 Summary 

Fig. 3 contains our Nf = 24 data and Nj = 17, 8, 4 and 2 data for N T = 4 - 8. 
Some of these data points correspond to real first order transitions, others describe just a 
crossover. For larger masses they correspond to the Z% transition, for smaller masses they 
describe the chiral transition. The agreement with the analytic prediction, especially with 
smaller Nj, is remarkable even for masses as small as m = 0.05 or below. The fact that 
the data appear to lie on a universal curve is a signal that the fermions induce an effective 
P whose strength is linear in Nj at fixed quark mass, down to very small mass. 

One can conclude that the effect of dynamical fermions for finite temperature transition 
is no more than an induced effective gauge coupling even for fairly small (m > 0.05) fermion 
masses. 
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5. Light and Heavy Flavors Together 



5.1 Nf = 2 + 1 

Several group studied QCD with two light and one heavy flavors. In these simulations 
the heavy fermion mass was between 0.1 and 0.25 and 4 to 20 times heavier than the light 
species. According to the previous chapter the effect of fermions with m = 0.1 — 0.25 on 
the finite temperature transition is well described by an effective action with induced gauge 
coupling given by Eqn. 6. Unless there is a strong interaction between the light and heavy 
flavors one would expect the same here. This hypothesis can be tested by comparing the 
shifts caused by a light species in the Nf = 2 and N j = 2 + 1 simulations assuming the 
effect of the heavy flavor can be described by Eqn. 6. Table VII shows this comparison 



(14) 



for recent simulations. Here we use the notation 

A(3f +1 = (/3? - (/3 r 2+1 + A/3^))/2 

A/3 2 = (/3?-/3 c 2 )/2. 

where (3$ is the Monte Carlo quenched critical coupling, /3 2 +1 is the Monte Carlo Nf = 2+1 
critical coupling and /3 2 is the Monte Carlo Nf = 2 critical coupling. Apanal j g 
analytically predicted induced gauge coupling due to the heavy fermions. If A/3 2+1 = A/3 2 , 
the shift in coupling due to the light quarks is independent of the presence of the heavy 
quark. If Eqn. 6 is valid for the light species, too, we expect A/3 2+1 = A/3 2 = A/3" nai The 
shifts from the Nf = 2 + 1 and Nf = 2 simulations agree within errors though A/3 2+1 is 
consistently smaller than A/3 2 which agrees with the analytic prediction Af3f nal . 

5.2 Nf = 2 + 8 

As another test of the interplay of light and heavy flavors, we performed simulations 
with two light flavors and eight heavy flavors on 8 3 x 4 lattices. 
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We were motivated to perform these studies by consideration of the deconfinement 
transition with Wilson fermions. Wilson fermions have 15 doublers for each light species. 
What is the role of the doublers? Do they only generate an effective gauge coupling or can 
they influence the low energy spectrum in a non-trivial way? Perturbatively, the 4 doublers 
sitting at the nearest edges of the Brillouin zone with one component of momentum equal to 
7r and three components of momentum equal to zero are much lighter than the others. It is 
plausible to assume that they give the most important contribution to the effective action, 
i.e. 2 flavors of Wilson fermions can be modeled as 2+8 flavors of staggered fermions. 
(Wilson fermions and 2+8 flavors of staggered fermions are of course not identical, since 
they have different flavor and chiral symmetry properties. This approach just models the 
effect of doublers.) 

For the heavy flavors we choose mass values = 0.88, 0.77, 0.665, and 0.4, corre- 
sponding roughly to the bare Wilson doubler masses atK = 0.17 — 0.21. The light masses 
were chosen as listed in Table VIII. 

These runs were performed on the Intel iPSC/860 hypercube at the San Diego Super- 
computer Center. The iPSC/860 and the code are described briefly elsewhere 7 . We used a 
truly hybrid algorithm for these simulations: the eight heavy flavors were simulated using 
the $ algorithm of Ref. 8, with a random noise term for the fermions which was refreshed 
at the start of each micro canonical trajectory. (The fermion fields were defined on all sites 
of the lattice to produce eight flavors.) The two light flavors were simulated using the R 
algorithm of Ref. 8; the noisy estimator for their determinant was updated throughout 
the simulation. We also performed a two flavor simulation at m = .04 for comparison. We 
used integration timesteps of At = 0.1 for the m; = 0.2 and 0.1 simulations. The smaller 
quark mass simulations were more sensitive to At systematics. We used At = 0.05 away 
from the transition for all the m = .04 simulations and switched to At = 0.02 near the 
transitions. 
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We display plots of the Polyakov loop and ipip for the light quark from our simulations 
with light quark mass 0.04 in Fig. 4. The heavy quark masses are 0.4, 0.665 and oo. The 
smaller step size points are shown as squares in the figure. The transition for the system 
with two light and eight heavy flavors appears to be much sharper than the transition for 
the system containing only two light flavors. It might be first order. Note that at this 
value of the light quark mass the Nf = 2 transition is a smooth crossover and the Nj = 8 
system does not have a first order transition for > 0.25 either. 

In the previous chapter we concluded that the gauge coupling induced by the heavy 
flavors is well described by the analytic formula. Using Eqn. 6 we compute the shift 
caused by one of the light flavors as in sect. 5.1 and compare it to the shift observed in 
the Nj = 2 simulations. We present these results in Table VIII. Our results for m/ = 0.1 
and 0.2 reproduce the Nj- = 2 simulation results and the analytic prediction, as we would 
expect following the successes recorded in the last chapter. Neither of our m; = 0.04 
results agree with the analytic formula. That could be explained simply as a breakdown 
of the analytic formula at light quark mass. However, with Nj = 2 or 4, simulations at 
rriq = 0.05 still agree with the analytic formula, as can be seen by comparing Tables V, 
VI, and the last entry of Table VII. What is even more surprising, the = 0.665 and 0.4 
data show a different shift in j3 from the same light quarks m; = 0.04. These facts, coupled 
with the qualitative sharpening of the transition at smaller m^, lead us to conclude that 
the eight heavy flavors have an observable influence on the light flavors in addition to an 
induced gauge coupling. The assumption that the heavy flavors are unimportant at low 
energies does not seem to hold. 

One might expect that this result would be even stronger if the light fermions are 
lighter. 

One might also expect similar behavior for Wilson fermions. In fact, one might expect 
an even stronger effect, since Wilson fermions include explicit interactions between the 
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light quarks and the doublers which are not present in this 2 + 8 flavor system. 
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6. /3 = limit 

128 flavors of fermions with mass m ~ 0.4 induce a gauge coupling 6/g 2 = (3^ ~ 7.6. 
That is large enough to deconfine an Nj< = 4 system even when the plaquette gauge 
coupling is zero. With large number of flavors one should see a confining-deconfining pure 
gauge phase transition in the (3 = limit as the function of the fermion mass. 

The naive analytical prediction in Sect. 2 predicted that for Nf > 16 flavors the 
fermions always decouple from the low lying gauge spectrum even in the (3 = strong 
coupling limit. Does that mean that for Nf > 16 at (3 = one will always find a deconfining 
phase transition for some value of the quark mass? 

We obviously cannot check this scenario numerically but we can study the iVy = 4 
finite temperature phase transitions in m at (3 = for different Nj values. Fig. 5a shows 
m cr it and 5b shows (3^ at the phase transition calculated from Eqn. 6 as the function 
of Nf. Since we observed strong metastability in all cases, we conclude that the phase 
transition with so many fermions is first order. 

The induced (3 lies in the range 7.5 to 8.5 for Nj > 80. The constancy of this result 
over a wide range of Nf indicates that the fermions do induce an effective gauge coupling 
which scales with Nj. This (3^ is not consistent with the quenched critical coupling 
(3c = 5.69 indicating that 6/g 2 : the coefficient of F^ V F^ V , does not equal (3 for small (3 
values. 

7. Conclusion 

We demonstrated that the effects of fermions on the finite temperature phase transition 
can be described by an induced effective plaquette term for masses as low as m ~ 0.05. 
The induced coupling is proportional to the flavor number and is independent of iVj>. The 
proportionality constant is given by a simple 1-loop formula. It is amazing that the simple 
formula for a fermion-induced shift in (3 works so well down to such small quark mass, 
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for degenerate mass fermions. From the point of view of lattice simulations of QCD, our 
results show that some dynamical quarks must be very light to cause interesting effects. 
A finite temperature simulation at some quark mass ought to show an induced (3 which 
is not given by the one-loop formula, before one could claim that a T = simulation at 
the same mass would be sensitive to the effects of dynamical quarks. This is just barely 
the case in contemporary dynamical fermion simulations. For example, the spectroscopy 
of the HEMCGC simulations at (5 = 5.6 with N f = 2 and m = 0.025 and 0.01 has 
been mapped onto quenched simulations at f3 = 5.935 and 5.95, respectively^ 4 These 
comparisons correspond to shifts per flavor of A/3 = .1675 and 0.175, respectively, to be 
contrasted with A(3 anal = 0.20 and 0.25 and finite temperature Monte Carlo shifts of about 
0.20 and 0.21. Thus they are in a regime where the sea quarks might be important for 
long distance dynamics. 

We simulated systems with 2 light and 8 heavy flavors to study the interaction of 
heavy and light quarks. For light masses m; > 0.1 we found no observable effect. For 
mi = 0.04 interaction with the heavy fermions as heavy as m = 0.665 can be observed in 
the finite temperature phase transition. We found the Nf = 8 + 2 transition is very sharp 
and its location cannot be predicted from the Nj = 2 transition assuming that the effect 
of the heavy flavors is described by an induced gauge coupling. These results may have 
applications to technicolor models, where one has to deal with the low energy effects of 
large numbers of heavy fermions as well as a small number of light fermions. Consequences 
of these results for Wilson fermions remain an open problem. 
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FIGURE CAPTIONS 



1. The real part of the Polyakov loop for m = 0.5, Nj = 24. 

2. Time evolution of the real part of the Polyakov loop for m = 0.5, Nf = 24. 

3. The induced gauge coupling divided by the number of flavors, A/3/Np from the sim- 
ulations described in this paper, compared with the curve from Eqn. 6, as a function 
of quark mass. Data are labeled with octagons for Nj = 24, pluses for Nj = 17, 
diamonds for Nf = 8, bursts for Nf = 4, and squares for Nf = 2. 

4. Plots of (a) the Polyakov loop and (b) ipip for simulations with two flavors of light 
quarks (m; = 0.04) and either nothing else (diamonds) or eight flavors of heavy quarks, 
of mass m/j = 0.665 (octagons and squares) or 0.4 (crosses and squares). The squares 
show data points from simulations with At = 0.02; all other data points used At = 0.05. 

5. (a) Plot of the location of the confinement-deconfinement transition at (3 = as a 
function of quark mass for several values of Nf. (b) The same data, but now interpreted 
as a function of the induced coupling inferred from the analytic expression. The line 
shows the location of the quenched deconfinement transition. 
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TABLE CAPTIONS 



I. A[3 as predicted by Eqn. 6 as the function of the quark mass. 

II. Nf = 24 simulations on 6 3 x 4 lattices performed by us. All the phase transitions, 
except the m = 1.00 one are very sharp, probably first order. At m = 1.00 there is 





only a broad crossover around (3 = 5.24. 


III. 


N f = 


17 simulations, from Ref. 9, at Nj< = 4. 


IV. 


N f = 


8 simulations, from Ref. 11 


V. 


N f = 


4 simulations performed by several groups 


VI. 


N f = 


2 simulations performed by several groups 


VII. 


N f = 


2 + 1 simulations from Ref. 2 and 1. 


VIII. 


N f = 


2 + 8 simulations performed by us. 
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Table I. A/3 as predicted by Eqn. 6 as the function of the quark mass. 



m 


.025 


.05 


.1 


.2 


.3 


.4 


.5 


.75 


1.0 


A/3 


.203 


.168 


.133 


.096 


.074 


.059 


.048 


.029 


.020 



Table II. N j = 24 simulations on 6 x 4 lattices performed by us. All the phase 
transitions, except the m = 1.00 one are very sharp, probably first order. At 
m = 1.00 there is only a broad crossover around /3 = 5.24. 



m 




A(3 MC /N f 


A(3 anal /N f 


1.00 


5.24(4) 


.0188(16) 


.0175 


0.75 


5.00(2) 


.0287(8) 


.0286 


0.60 


4.76(2) 


.0387(8) 


.0388 


0.50 


4.62(2) 


.0446(8) 


.0478 


0.25 


3.90(5) 


.0746(20) 


.0840 



Table III. N f = 17 simulations, from Ref. 9, at N T = 4. 



m 


Pc 


Ap MC /N f 


A(3 anal /N f 


0.50 


5.025 


.039 


.0478 


0.25 


4.6(1) 


.064 


.084 


0.10 


4.3(1) 


.082 


.133 



Table IV. Nf = 8 simulations, from Ref. 11 



m 


Pc 


A/3 MC /N f 


A(3 anal /N f 


N T 


1.0 


5.54 


.0187 


.0175 


4 


0.50 


5.31 


.0475 


.0478 


4 


0.25 


5.025 


.083 


.084 


4 


0.10 


4.80(1) 


.111(1) 


.133 


4 


0.10 


4.95 


.12(1) 


.133 


6 


0.05 


4.75 


■14(1) 


.168 


6 
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Table V. Nf = 4 simulations performed by several groups. 



m 


Pc 


Af3 MC /N 


A[3 anal /N f 


± 


Ref. 


0.5 


5.50 


.04 


.048 


4 


13 


0.5 


5.45 


.055 


.048 


4 


15 


0.4 


5.42 


.063 


.059 


4 


15 


0.3 


r or 

5.35 


.08 


.074 


4 


13 


0.2 


5.255(5) 


.104(2) 


r\(\ r o 

.0958 


4 


16 


0.1 


5.130(5) 


.136(2) 


.1334 


4 


16 


0.05 


5.04 


.163 


.168 


4 


17 


0.0375 


4.99 


.173 


.183 


4 


18 


0.0375 


5.02 


.168 


.183 


4 


17 


0.025 


4.98(2) 


.175 


.203 


4 


16 


0.0125 


4.919 


.190 


.240 


4 


18 


0.01 


4.95 


.185 


.25 


4 


17 


0.25 


5.509 


.090 


.084 


6 


19 


0.10 


5.322 


.137 


.133 


6 


19 


0.075 


5.25 


.155 


.149 


6 


19 


0.065 


5.22 


.162 


.155 


6 


19 


0.05 


5.175 


.174 


.168 


6 


19 


0.025 


5.130 


.187 


.203 


6 


17 


0.01 


5.08 


.199 


.25 


6 


17 


0.025 


5.25 


.188 


.203 


8 


20 


0.01 


5.15(5) 


.213 


.25 


8 


21 
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Table VI. Nf = 2 simulations performed by several groups. 



777 
III 


/?- 
yc 




Aoanal / at- 

L-i/J I f 


Nrr 


J. LCI. 


1.0 


5.63 


.02 


.018 


4 


13 


0.4 


5.54 


.065 


.059 


4 


13 


0.2 


5.48 


.095 


.096 


4 


13 


0.1 


5.38 


.15 


.133 


4 


15 


0.05 


5.34 


.165 


.168 


4 


13 


0.025 


5.2875 


.197 


.203 


4 


15 


0.0125 


5.271 


.21 


.24 


4 


22 


0.01 


5.265(10) 


.212 


.25 


4 


2 


0.025 


5.445 


.212 


.203 


6 


23 


0.0125 


5.42(1) 


.225(10) 


.239 


6 


23 


0.0125 


5.5375 


.23 


.239 


8 


24 
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Table VII. iVy = 2 + 1 simulations from Ref. 2 and 1. 



mi 


™h 






A/3 2+1 


A/f 






0.025 


0.025 


5.132(2) 


.20 


.18 


.20 


.20 


4 


0.025 


0.10 


5.171 


.13 


.20 


.20 


.20 


4 


0.0125 


0.25 


5.199(2) 


.084 


.20 


.23 


.24 


4 


0.00833 


0.1667 


5.325(25) 


.11 


.22 




.26 


6 



Table VIII. Nj = 2 + 8 simulations performed by us. 



mi 


™>h 




^panal 


A(3f +8 






0.04 


0.665 


5.065(5) 


.271 


.170(3) 


.183(5) 


.18 


0.04 


0.40 


4.89(1) 


.474 


.155(5) 


.183(5) 


.18 


0.10 


0.77 


5.18(1) 


.22 


.138(5) 


.15 


.133 


0.20 


0.88 


5.275(25) 


.18 


.11(2) 


.095 


.096 


0.04 


oo 


5.275(25) 







.183(5) 


.18 
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